Proximity induced pseudogap in mesoscopic superconductor/normal-metal bilayers 
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Recent scanning tunneling microscopy measurements of the proximity effect in Au/La2-a;Sr x Cu04 
and Lai.55Sro.45Cu04/La2- a: Sr a; Cu04 bilayers showed a proximity-induced pseudogap [Yuli et 
al., Phys. Rev. Lett. 103, f 97003 (2009)]. We describe the proximity effect in mesoscopic 
superconductor/normal-metal bilayers by using the Bogoliubov-de Gennes equations for a tight- 
binding Hamiltonian with competing antiferromagnetic and d-wave superconductivity orders . The 
temperature dependent local density of states is calculated as a function of the distance from the 
interface. Bound state due to both d-wave and spin density wave gaps are formed in the normal 
metal for energies less than the respective gaps. If there is a mismatch between the Fermi velocities 
in the two layers we observe that these states will shift in energy when spin density wave order is 
present, thus inducing a minigap at finite energy. We conclude that the STM measurement in the 
proximity structures is able to distinguish between the two scenarios proposed for the pseudogap 
(competing or precursor to superconductivity). 
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PACS numbers: 74.45+x, 74.50. +r 



One of the most important issues in cuprate high tem- 
perature superconductors is the existence of a pseudo- 
gap (PG) phase in the underdoped region 1,2 . It has re- 
ceived heavy debate about its nature as its origin may 
hold key information regarding the pairing mechanism 
in these materials. One paradigm assumes that the PG 
is only a manifestation of a competing or coexisting or- 
der of superconductivity and it has no direct relation to 
the pairing^—. However, other scenarios regard the PG 
as a superconducting precursor stated—, with close re- 
lation to Cooper pairing. Despite many studies on the 
subject, it remains unclear whether the PG phenomenon 
is related to superconductivity or not. 

Very recently, the temperature evolution of the 
proximity-induced gap in Au/La2_ 3; Sr 2 ;Cu04 and 
Lai.55Sro.45Cu04/La2-a;Sr x Cu04 bilayers was reported 
by Yuli et al. 12 using scanning tunneling microscopy 
(STM) and spectroscopy (STS) measurement to exam- 
ine whether unique spectral properties associated with 
pairing are present above the superconducting transition 
temperature T c in the PG temperature regime. For 
bilayers comprising an underdoped La2- a; Sr : j : Cu04 
(LSCO) film, there exists a smooth evolution of the 
superconductor proximity gap into a proximity-induced 
PG, as the temperature was raised above T c . In contrast, 
the proximity gap disappeared close to T c in bilayers 
comprising an overdoped LSCO layer. They claimed that 
the similar spatial dependence of the proximity induced 
PG and the superconductor proximity gap indicates 
that the origin of the PG is related to superconductiv- 
ity. How to understand these behaviors theoretically 
is still a challenging question. The proximity effect 
between a superconductor and a normal metal has 
been thoroughly studied using various techniques^—. 
However, an investigation of the proximity effect in 
normal-metal/superconductor bilayers above T c is still 
lacking. 

In the present paper we use numerical solutions of 



the Bogoliubov-de Gennes (BdG) equations based on 
a model Hamiltonian with competing d-wave supercon- 
ducting (DSC) and antiferromagnetic (AFM) orders at 
finite temperature and examine the proximity-induced 
gap in mesoscopic normal-metal/superconductor bilay- 
ers. In cuprate superconductors, the spin density wave 
(SDW) order with stripe modulation and the accompa- 
nying charge order may emerge at high temperature T, 
and could exist above the DSC transition temperature 
Tj£. Our numerical analysis focuses on the interplay 
between DSC and SDW orders when in proximity with 
the metal as well as the local density of states (LDOS) 
as a function of temperature and distance from the inter- 
face. States induced in the normal metal region depend 
on both DSC and SDW orders present in the supercon- 
ducting side and on the Fermi velocity mismatch at the 
interface. Our results provide important new insight into 
the interpretation of recent STM experiments 12 . 

In order to describe the superconductor/normal-metal 
bilayers we use the tight-binding extended Hubbard 
Hamiltonian by assuming that the on-site repulsion Ui 
is responsible for the competing AFM order and the 
nearest-neighbor attraction Vy for the DSC pairing in 
the superconducting region: 
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where ty = t are the nearest-neighbor hopping integral. 
c\ a (c[ CT ) are destruction (creation) operators for electron 

of spin <t, n\ a = c \ a c -\a is the number operator, and [i is 
the chemical potential determining the averaged electron 
density n. The SDW and DSC orders have the following 
definitions respectively: A? DW = Ui{cLa-\ — cjiCy.) and 
Aij 
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formation, 



- CL|Cjf)/2. Using the Bogoliubov trans- 
En [Uialna-w?* jig], the Hamiltonian 





FIG. 1. The spatial variation of DSC order parameter 
Af (a) and staggered magnetization Mf (b) as a function 
of N x for different temperatures T in a mesoscopic normal- 
metal/superconductor bilayer. "S" and "N" denote respec- 
tively the superconducting region and normal metal. The 
parameter values are: U — 2.5, ft = 0.9, and t — 0. 



in Eq. (1) can be diagonalized by solving the resulting 
BdG equations self-consistently, 
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where Hij a = — iy 



[Ui(nia) — A*]^y • Using open bound- 
ary conditions we obtain the eigenvalues {E n } with eigen- 
vectors {u",v™}. 



The self-consistent conditions are: 
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where /(.E) = (e E l ksT + 1) x is the Fermi-Dirac distri- 
bution function. The DSC order parameter is defined at 
site i as Af = (A i+ex .i + Aj_ ex .j - A u+ey - A u _ ey )/4, 
where e X) ( y ) denotes the unit vector along the x, (y) direc- 
tion. Throughout this work, the distance is measured in 
units of the lattice constant a, and the energy is scaled to 
4. In the numerical calculations, we take ks = a = t = 1 
for simplicity. The averaged electron density is chosen as 
n=0.9, corresponding to the underdoped level (5=0.1. For 
an appropriate initial order parameter profile, the Hamil- 
tonian is numerically diagonalized and the obtained elec- 
tron wave functions are used to calculate the new param- 
eters for the next iteration step. The solution is found 



FIG. 2. The spatial variation of Af (a) and M° (b) as a 
function of N x for different temperatures T in a mesoscopic 
normal-metal/superconductor bilayers. The parameter values 
are: U = 2.5, ft = 0.9, and t = -0.2. 



when the relative error in the gap function between suc- 
cessive iterations is less than the desired accuracy. 



The LDOS of the energy E at the position i can be 
written as 

2JV 

Pi (E) = Y,M\ 2 $(E-E n ) + \v?\ 2 6(E + E n )], (6) 



where 8(E) is the Dirac delta-function, which is broad- 
ened into —- 



- (E-E') 2 +r 2 ' w ith r = 0.01 being chosen. 
Pi(E) is proportional to the local differential tunneling 
conductance which could be measured by STM experi- 
ments. 

First, we consider mesoscopic superconductor / 
normal-metal bilayers with sizes of iV a; xiVj / =32x32 
(12x32) for the superconducting (normal) region. The 
spatial variation of the DSC order parameter Af (a) and 
the staggered magnetization Mf (b) along the x axis for 
various temperatures T are plotted in Fig. Q] for chosen 
parameter values: U = 2.5t and V = 1.24 and 4' = 0. 
Here, the staggered magnetization of the induced AFM 
or SDW order is defined as M? = (-l) l (n^ -n i]r ). Fig. M 
shows the DSC and SDW order parameters for the same 
pairing potentials but for a different next-nearest neigh- 
bor hopping 4' = —0.24 in order to properly account for 
the band structure of high-T c cuprates. One can easily 
see that, for both cases, the y-oriented stripe modula- 
tions of DSC and coexisting SDW orders are present at 
lower temperatures. Note that, the oscillation periods de- 
crease with increasing the next-nearest neighbor hopping 
strength, while the period of the SDW is almost half of 
that of the DSC. Moreover, the periodicity of the strip- 
modulated DSC and SDW orders for 4' = seem not 
to be sensitive to T. When the temperature increases, 
the DSC order is suppressed and finally disappears at 



T c . Further increasing T, the staggered magnetization 
decreases to zero at the Neel temperature Tjv =0.3. In 
this case the nature of the pseudogap is the competing 
SDW order present only at low doping and higher tem- 
perature. This coexistence could be the reason for the 
modulation of the superconducting gap observed in STM 
experiments on underdoped cuprates 19 . 

Shown in Figs.Q]and[l both DSC (panel (a)) and SDW 
(panel (b)) leak into the normal region due to the prox- 
imity effect. The usual exponential decay is observed, 
with a leaking distance which is of the order of the lat- 
tice constant. Similar to bilayers made of conventional 
SC and normal metals, bound states will be induced in 
the normal region for energies less than the supercon- 
ducting gap. In order to see this we plot in Fig. [3] the 
LDOS around the SC/normal metal interface. Similar 
to Ref. [l8| we find that in the superconducting region as 
the temperature is increased above T c a gap still exists 
due to the coexisting SDW order. Unfortunately, due 
to the discreteness of the energy spectrum and the finite 
temperature smearing, we cannot draw meaningful con- 
clusions from the LDOS calculated in the normal-metal 
region. The main observation from Fig. [3] is the modifica- 
tion of the LDOS in the normal metal region for energies 
below the DSC or SDW gaps. 

In order to better understand the formation of bounds 
states in the normal region we performed a slightly dif- 
ferent calculation. With the insight given by the self- 
consistent calculation done on a smaller sized system, we 
calculated the LDOS for a bilayer of size N x x N y = 
300 x 300, with a normal metal region of size L = 4. 
While still using the same model Hamiltonian we do 
not recalculate the order parameters self-consistently but 
consider them as step functions. If we are interested only 
in the LDOS at the surface of the normal-metal this ap- 
proximation is reasonable since all the order parameters 
are set to zero in the normal region while if the coher- 
ence length is small it results in a sharp drop of the order 
parameters at the SC/N interface. Now instead of diago- 
nalizing the full Hamiltonian we approximate the Green's 
functions and hence the LDOS in terms of Chebyshev 
polynomials. The procedure is presented elsewhere 2 ^. 

Similar to an Andreev reflection at the N/S interface 
for quasi-particles with energies lower than the supercon- 
ducting gap, at the N/AFM interface a peculiar reflec- 
tion (Q-reflectio n 21 i 22 ) occurs if the energy is less than 
the SDW gap. Quasi-particles with momentum hp are 
scattered to states with momentum k^p + Q where Q is 
defined by the AFM order. If there is a potential barrier 
at the interface or a mismatch in the Fermi velocity, in 
addition to the Q-reflection, specular reflection also oc- 
curs thus influencing the formation of the bound states. 

The main point in the STM measurement^ was the 
fact that the observed gap is always at zero energy 
whether the normal metal layer is overdoped LSCO or 
Au. To simulate these conditions we consider three situ- 
ations: no Fermi velocity mismatch, Fermi velocity mis- 
match from a shift in the chemical potential in the nor- 
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FIG. 3. The LDOS for positions around the interface for 
t = -0.2* with T = 0.1, T c = 0.17 and T = 0.2. The other 
parameter values are: U — 2.5 and n = 0.9. The thick black 
line represents the LDOS at the interface while lower(upper) 
lines are in the SC(normal) region. The normal state LDOS 
is subtracted for all lines in order to suppress oscillations due 
to the interface with vacuum. 



mal region (5fi = 0.5) or Fermi velocity mismatch from 
a different band structure (*' = in the SC region while 
t' = —0.2t in the normal region). In Fig. 2] we plot the 
calculated LDOS for three choices of the order parame- 
ters. First, in Fig. gja) only DSC (A d = 0.4*) is consid- 
ered. Subgap structures appear for all three choices of 
interfaces but the minigap is independent of this choice. 
Second, in Fig. g^b) only SDW {A SD w = 0.4*) is con- 
sidered. Subgap features are also observed but depend- 
ing whether there is Fermi velocity mismatch or not, the 
bound states shift in energy. The minigap now shifts to 
higher energy. Third, in Fig. SKc) we consider a coex- 
istence of DSC and SDW (A d = 0.4*, A SD w = O.lt). 
In this case a minigap is still observed when the next- 
nearest neighbor hopping is changed in the normal layer 
although the particle-hole symmetry is broken. When 
the chemical potential is shifted in the normal metal, a 
suppression of the LDOS is still observed but the gap is 
closed. This behavior is not consistent to the experimen- 
tal STM findings^ 2 -, which showed that the minimum of 
the proximity induced pseudogap is always at the Fermi 
level irrespective of the normal metal type (overdoped 
LSCO or Au). In order to explain the observed prox- 
imity induced pseudogap only scenarios which involve 
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FIG. 4. The LDOS at the surface of the normal metal region 
for three choices of order parameters: (a) Ad = 0.4£, Asdw = 
0.0, (b) A d = 0.0, Asdw = OAt, (c) A d = OAt, A S dw = 
O.lt. The lines represent three Fermi velocity mismatches at 
the N/SC interface: no mismatch, mismatch due to change 
in chemical potential and mismatch due to change in next- 
nearest neighbor hopping. 



reflections with Q = 0, e.g. fluctuating SC order pa- 
rameters or loop current orders^, should be considered. 



In conclusion, we have investigated the proximity ef- 
fect in mesoscopic normal-metal/superconductor bilayers 
by numerically solving the BdG equations based on an 
effective model Hamiltonian with competing AFM and d- 
wave superconductivity interactions. The self-consistent 
solution gives oscillating DSC and SDW order param- 
eters with the period of oscillations determined by the 
value of the next-nearest neighbor hopping amplitude. 
For low doping there exists a temperature range above 
T c in which only the SDW order is stable. As expected 
we observe in the metallic region the proximity induced 
pair correlations and the localized bound states for ener- 
gies less that the DSC (or SDW) gap. More importantly 
by using a larger system we were able to confirm the con- 
clusion presented in Ref. [12| which argues that an SDW 
order above T c would shift the minimum gap energy if the 
normal metal is changed from overdoped LSCO to Au. 
By using a model microscopic Hamiltonian we showed 
that a mismatch in the Fermi velocity at the N/SC in- 
terface would shift the induced minigap. Depending on 
the nature of this mismatch (potential barrier or different 
band structure) the shift should also be observed below 
T c in the temperature range where DSC and SDW order 
coexist. 
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